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The widely used Gross-Pitaevskii equation treats only co- 
herent aspects of the evolution of a Bose-Einstein condensate. 
However, inevitably some atoms scatter out of the condensate. 
We have developed a method, based on the field theory formu- 
lation, describing the dynamics of incoherent processes which 
are due to elastic collisions. We can therefore treat processes 
of spontaneous emission of atoms into the empty modes, as 
opposed to stimulated processes, which require non-zero ini- 
tial occupation. 

In this article we study two counter-propagating plane 
waves of atoms, calculating the full dynamics of mode oc- 
cupation, as well as the statistics of scattered atoms. The 
more realistic case of Gaussian wavepackets is also analyzed. 

In recent years, there is a growing number of exper- 
iments in which the atomic Bose-Einstein condensate 
undergoes a nontrivial dynamics. Some examples are 
bouncing condensates [Q, dark solitons in a condensate 
[g,||, vortices and their dynamics Q|, Bragg splitting of 
the condensate Q and many others. 

A remarkably universal tool describing all these exper- 
iments is the celebrated Gross-Pitaevskii equation (GPE) 
[M,R|. It describes a coherent evolution of the mean 
atomic field. In the most common interpretation, its time 
dependent version assumes that each atom of the system 
undergoes identical evolution. This is a good assump- 
tion since in these experiments, initially the wave-packet 
of the system contains thousands of particles. To use 
a term borrowed from quantum optics, the time depen- 
dent GPE describes stimulated processes. At first glance, 
atomic four wave mixing looks like an exception. It is 
not, since in this case the nonlinear term in the GPE 
contains a macroscopic source term of the fourth wave. 
The analogous optical process is also described well by 
the c-number version of Maxwell equations with a non- 
linear polarization ||9|. Similarly, the second harmonic 
generation requires no field quantization, but the reverse 
process of the parametric down conversion starts-off due 
to vacuum fluctuations and requires quantized photons 
P|. Also, the mechanism of spontaneous emission is nec- 
essary for single atom radiative decay and its multi-atom 
generalization known as superfluorescence |Q . 

In the process of a collision between two condensates, 
inevitably, some atoms from both droplets of quantum 
matter would scatter away from the two condensates. 
This is a loss mechanism which is not accounted for by the 
conventional GPE. The models were developed in which 
the effects of this loss on the condensates dynamics were 



incorporated in the GPE within the slowly varying en- 
velope approximation technique [Q. It is a purpose of 
this article to propose a simple method of calculating not 
only these losses but also quantum statistical properties 
of the atoms scattered away. Our starting point will be a 
second quantized version of the many-atom Hamiltonian. 
Taking our inspiration from quantum optics, we shall ap- 
proximate this Hamiltonian in such a way, that the only 
quntized modes of the atom field will be those which ini- 
tially are in the vacuum state. Solving the Heisenberg 
and Schrodinger equations for the empty modes we shall 
be able to answer all relevant questions concerning the 
lost atoms. 

The Hamiltonian describing a multi-atom system mov- 
ing in a large box of length L and undergoing cold colli- 
sions via a contact potential reads 

H= d^r^\Y,t)—^{Y,t) + 

Jv=L^ 2m 

+ f/ d''r¥{v,t)¥{r,t)^ir,t)^{v,t) (1) 

where ^ denotes an atomic field operator satisfying the 
standard equal time Bose commutation relations and 
coupling constant g = Anh a/m, where a is the s-wave 
scattering length. 

In this article we consider a collision processes of two 
large Bose-Einstein condensates. In particular we want 
to compute the part of the atomic field consisting of 
atoms scattered away from both condensates. In this 
process two out of four field operators in the last term 
in (|l|) may be replaced by the c-number colliding wave- 
packets. In the scattering process two atoms from these 
wave-packets are destroyed and two atoms appear in the 
initially empty modes. As the losses are assumed small, 
in the crudest approximation we do not take into account 
the modification of the colliding condensates due to the 
scattering process. In this way we obtain a Hamiltonian, 
which for the problem considered here resembles the ra- 
diation of given current in quantum electrodynamics |12| ] 
or a parametric down conversion with undepleted pump 
in quantum optics ||ll[ . Its hermitian form contains both 
creation and annihilation terms: 
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H= I d'r^lir,t)j-^,{r,t) 



d^r^lir, i)*l(r, t)^bQir, t)V'-Q(r, t) + h.c. (2) 



Here V'±Q are the wave functions of colliding conden- 
sates (in the center of mass reference frame) and ^e is 



the field operator of the empty modes. Hamiltonian (g) is 
quadratic in the atomic field operators. Hence, the field 
equations will be linear and tractable. We pay a price 
for this dramatic simplification: unlike (0) the Hamilto- 
nian (0) does not conserve the number of particles. Also 
the dynamics of the condansates cannot be found in this 
approach, GPE must be solved beforehand. 

The roots of this approximation, together with the 
validity conditions, might be traced back to the time- 
dependent generalization of the Bogoliubov approach, in 
which the boson field operator is decomposed into the 
macroscopically occupied part (i.e. the condensate) and 
the quantum corrections. As our interest lies in collision 
processes of Bose- Einstein condensates, the wavefunction 
of the whole system is just the sum of colliding conden- 
sate's wavefunctions. The Bogoliubov decomposition can 
be therefore written as: 



*(r, t) = VQ(r, t) + Vj_Q(r, t) + #e(r, t) 



(3) 



Inserting (g) into the Hamiltonian (|l|) one obtains a 
collection of different terms, of which only two remain 
in Hamiltonian (^. We shall argue that such an ap- 
proximation gives correct results if and only if the ki- 
netic energy associated with the center of mass motion 
is much larger than the interaction energy per particle, 
h Q^ /{2m) ^ gn, where n is the average density of the 
particles in the condensates. This condition is readily ful- 
filled in current experiments (see for example |l3-n3]). 

We can prove the above statement in the simplest case 
of two counter-propagating plane matter waves. The 
wavefunctions of the condensates evolving according to 
the free Hamiltonian are the following: 



V'±Q(r,i) 



V^^e^'Q'-e- 



ihi 



(4) 



where n±Q is the density of particles. Let us assume for 
simplicity that both waves are equally populated, nq = 
n_Q = n. In a box of length L, the boson field operator 
can be decomposed into normalized plane waves: 



vI/e(r,t)-L-3/2^e-<i>-aq(i) 



(5) 



where aq(t) is an annihilation operator of an atom of 
the wavevector q. Bogoliubov approximation implies the 
following: 
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This way a set of coupled linear equations for cre- 
ation/annihilation operators is obtained. To solve it we 
have introduced a cut-off in the number of modes taken 
into account and then solved the resulting finite system 



algebraically. The numerically found populations of dif- 
ferent modes are plotted in Fig. |l|; the figures differ by 
the value of total kinetic energy of colliding waves. 
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FIG. 1. Bogoliubov approximation: population of modes 
to the power 1/3 (to make secondary rings visible) on a mo- 
mentum plane. Two initial plane waves are located at (Q, 0) 
and (— Q, 0). The parameter r is the ratio of the kinetic energy 
associated with the center of mass motion to the interaction 



energy per particle r 
our approximation. 
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Large r are well described by 



Different terms in Eq. (p) result in different aspects of 
mode's population in Fig. U The terms proportional to 
—ih-^ap and alp (the only ones that remain in our ap- 
proximation presented in this paper) are responsible for 
the elastic scattering of atoms from two opposite waves. 
They form the main circle centered at p = in nromen- 
tum plane that appears on the graphs. Increasing the 
value of kinetic energy (increasing r) results in relative 
sharpening of the circle p k, Q. One can see that for 
large values of Q populating this circle is the main effect 
of spontaneous emission. 

The other terms, — 2i^ (ap-2Q + 2ap + ap+2Q) are 
responsible for the frequency shifts, i.e. the fact, that 
the circles' radii are smaller than Q. Due to these terms 
also anisotrophy in the circle population appears (scat- 
tering perpendicular to the axis of the colliding waves 
is the most probable). However, these effects are get- 



ting negligible as Q increases - the circles grow, finally 
reaching p = Q and flatten. 



The last terms, -i^e^^'^'^'*/™ (a2Q-p 
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are responsible for secondary rings, which are a result 
of processes in which two atoms from the same matter 
wave scatter. After the collision one of them is placed 
in the main circle (because of the Bose enhancement) 
and the other must go to the secondary ring to fulfill the 
momentum conservation. However, the kinetic energy in 
such a process is not conserved: the energy of atoms after 
the collision is greater than before; the difference is taken 
from the interaction energy. Thus, as the ratio of kinetic 
energy to the interaction energy increases (parameter r) 
these processes are less and less important and finally 
practically disappear for large Q. 

Hence we have shown that increasing the mutual veloc- 
ity of the condensates results in sharpening the resonance 
condition and substantial depopulation of the modes with 
p^Q. 

Let us now consider the situation of two counter- 
propagating plane matter waves within our approxima- 
tion in more detail. Substituting (0) and ^ into Hamil- 
tonian (0) we obtain: 



^ - E (l^^w^'^w- 



^^e--^*<(t)at^(t) + ..c) (7) 



Clearly, the Hamiltonian of the system conserves momen- 
tum: two atoms from colliding plane waves of momenta 
hQ and —hQ scatter into atoms of momenta Tiq and 
— ?iq, the total momentum remains zero. 

The Heisenberg equations of motion for annihila- 
tion/creation operators can be solved analytically. The 
solution takes the following form: 

apit) = e-'^^* (^i(i,p)ap(0) - F2(t,p)alp(0)) (8) 
where: 
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and A(p) = (^^^) - i^^^^) ■ Modes can be 
divided into two classes according to their evolution. The 
ones which satisfy the modified energy conservation law, 
i.e. the kinetic energy difference between the incident and 
scattered atoms is smaller then the mean field energy: 
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(10) 



have their occupation growing exponentially with time, 
which is a Bose enhanced process. This relaxation of the 
energy conservation has an analog in quantum optics, 
known as power broadening Pq]. The modes that do 
not satisfy the resonance condition dTo) oscillate with 
the population of order unity. 

It follows from (ph that during the dynamical evolution 
the state of the system remains a tensor product of states 
from 2-mode subspaces. More specifically, if at i = only 
modes Q and — Q were populated, each by N/2 atoms, 
then at subsequent times the state of the system is of the 
form: 



|7V/2,7V/2)q® ( ^mt))^ 
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(12) 



and |fc,A:)q denotes a state with k atoms of momentum 
h(\ and k atoms of momentum —%(\. Hence, though the 
state of the system remains pure during dynamics, the 
statistics of the scattered atoms is geometric (the same 
as in a thermal state, which is mixed). Of course our 
solution is just an approximation to the particle conserv- 
ing one, in which the scattered atoms would be entangled 
with the highly occupied modes with suitable number of 
particles missing. 

To gain deeper insight into statistical properties of 
scattered atoms we calculate the ff^^HPil) correlation 
function: 
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(13) 



Again under the assumption that at t = all modes ex- 
cept the ±Q are empty the correlation function becomes: 



5^'Hp,q) = i 
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(14) 



We have atom bunching (g*-^' > 1) for two cases: p = — q, 
which is a result of the momentum conservation law and 
p = q, which refiects the Bose enhancement. Thus, it is 
expected that during the collision, scattered atoms will 
come out in the form of pairs of spikes, randomly placed 
at the p = Q circle. Only after the averaging over many 
shots the isotropy will be recovered. 

To get closer to reality, let us analyze the situation 
of two atomic wavepackets colliding. At t = the 
wavepackets envelopes are /±Q(r) and move with mo- 
menta ±?iQ. Then, at time t the wavefunction of each of 
the packets takes the form: 



V^±Q(r,t)=L--''^/±Q(k±Q)e- 



ikr — iS,^ 



(15) 



where /±Q(k) denotes the Fourier transform of /±Q(r). 
We are using wavef unctions rather than a density ma- 
trix since we are deahng with macroscipically populated 
matter waves. And even further simpUfication is made 
by replacing GPE solution by a free evolution of both col- 
liding wavepackets. We have therefore entirely neglected 
the mean field energy in comparison with the kinetic en- 
ergy of the center of mass motion. 

Hence, the hamiltionian (0) can be written as: 
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(16) 



In general the resulting Heisenberg equations for annihi- 
lation/creation operators cannot be solved analytically, 
but can be treated perturbatively The first order in the 
coupling constant g gives the following result: 
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To evaluate total spontaneous emission losses, one has to 
consider the quantity: 



5(t)=$^(at,(t)ap(t)) 



(18) 



where the average is taken over the vacuum. 

To calculate ( p^ ) we will make further specifications. 
We will consider Bragg splitting of a condensate (some- 
times called half-collision), which is easier to realize in 
experiment than a full collision between two condensates 
Q . We assume that the initial wavepackets are the Gaus- 
sians: 



/q(i") = /-Q(r) =Kexp 
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(19) 



where k fulfills the normalization condition (the colli- 
sion is along the x-axis). We will also assume non- 
dispersive evolution of wavepackets by neglecting the 
term exp(i?iuvT/m) in (|l^), which is justified by the 
fact that the scattering takes place while the wavepackets 
overlap, and for this short time scale the dispersion does 
not play a substantial role. In other words, characteristic 
time scale of the scattering process is assumed to be very 
short in comparison with both nonlinear and dispersive 
time scales. 

Within these approximations the total spontaneous 
emission losses ( |l^ can be evaluated analytically: 



S{t) = W^ -^ ^^ / (fp / dr / dr' 
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In fact, one more integral can be performed analytically, 
but the resulting formula is even more complex. It is also 



worth noting that the order of integrations in (20) must 
be kept - otherwise one can easily arrive at a divergent 
expression. 

In the case of spherical symmetry of the colliding 
wavepackets {ax,y,z ='-'^) the angular distribution of scat- 
tered atoms is isotropic and the integral with respect to 
the angles of p can be performed. We have noticed that 
effectively the total losses do not depend on |Q|, provided 
that the momentum is sufficiently large to produce inter- 
ference fringes in the splitting wavepackets: Qa 3> 1. 
We have checked that in this range total losses are pro- 
portional to: NiN2{a/a)'^. For numerical calculations 
we assumed A^i = 7V2 = 5 • 10^ sodium atoms in each 
wavepacket of width a = 30^m and the result is plotted 
in Figure |[ For these parameters the value of total spon- 
taneous emission losses is below 1% of the total number of 
particles. The curve is universal if plotted against dimen- 
sionless time t = t hQ/{mcr). In this graph we also show 
for comparison elastic scattering loss from the condensate 
calculated from GPE with complex scattering length [0 , 
where we neglected dispersion and depletion of the col- 
liding wavepackets. The latter may be described by an 
analytic formula: S{t) = 4iVi7V2(a/cr)2Erf(t\/2) (shown 
as a dashed line in Figure 0, but the curves are practi- 
cally indistinguishable). Both methods are first order in 
scattering length and agree extremely well. 
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FIG. 2. The growth of spontaneous emission losses relative 
to the total number of atoms versus time in units ma/{hQ). 



The influence of the geometry of the colliding 



wavepackets on the total losses has been investigated for 
cylindrically symmetric gaussians. For a wide range of 
the aspect ratio the total losses do not depend on the 
geometry (provided that the condition Qax ^ 1 is ful- 
filled). Similar effect is predicted by the complex scat- 
tering length method. 

The method presented here shows extraordinary good 
agreement with the predictions of the complex scatter- 
ing length technique within the approximations discussed 
earlier. It can be shown (using dynamic equations for non 
commutative cumulants) that within first order approx- 
imation both methods are essentailly equivalent [jl^. It 
would be interesting to go beyond these approximations 
to check the predictions of [^ , we are currently studying 
the Bose enhancement for realistic condensate droplets. 
It is worth stressing however that our method not only 
measures total spontaneous emission losses during the 
collison, but also gives insight into the statistical proper- 
ties of the scattered atoms. 

The recent experiment by MIT group [[l5| , measuring 
the losses in the splitted condensate, was performed in 
different regime than discussed in this article. The mea- 
sured fraction of collided atoms was up to 20%, which 
is clearly beyond the first order perturbation theory de- 
scribed above. Authors of ||l^] argue that the increase of 
the fraction of collided atoms with the number of out- 
coupled impurities (in our language that is the total loss 
versus the number of atoms in one of the packets, keep- 
ing the total number of atoms in both condensates fixed) 
is an evidence of collective self-amplification and is due 
to bosonic enhancement. It is striking however that the 
shape of this dependence (Fig. 3 of |l5|) resembles the 
parabola obtained from the formula (pfl) . 

To summarize: We have proposed a method that al- 
lows to calculate the spontaneous emission losses of col- 
liding condensates. Contrary to earlier attempts [Q, our 
approach uses wave functions rather than particle den- 
sities. The coherence of the colliding droplets is there- 
fore taken into account. The method is applicable to 
sufficiently quickly moving condensates, i.e. when the 
spread of momenta in each condensate is much smaller 
than their relative momentum and in this regime is essen- 
tially equivalent to Bogoliubov approach. For the sim- 
plest case of colliding plane waves we pointed out the 
thermal-like multiplicity distribution of scattered atoms 
in each mode and the bunching reflecting the momentum 
conservation on one hand and the bosonic enhancement 
on the other. Similar effects must be also present for the 
colliding wavepackets. 
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